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CS] ■ ABSTRACT 

I I , Context. In the past decade, much progress has been made in characterising the processes leading to the enhanced deuterium fractionation observed 

■^L in the ISM and in particular in the cold, dense parts of star forming regions such as protostellar envelopes. Very high molecular D/H ratios have 

r H\ been found for saturated molecules and ions. However, little is known about the deuterium fractionation in radicals, even though simple radicals 

^*, often represent an intermediate stage in the formation of more complex, saturated molecules. The imidogen radical NH is such an intermediate 

r^ species for the ammonia synthesis in the gas phase. Many of these light molecules however have their fundamental transitions in the submillimetre 

O I domain and their detection is hampered by the opacity of the atmosphere at these wavelengths. HerschelfHWl represents a unique opportunity to 

I study the deuteration and formation mechanisms of species not observable from the ground. 

O Aims. We searched here for the deuterated radical ND in order to determine the deuterium fractionation of imidogen and constrain the deuteration 
,y mechanism of this species. 

C/2 Methods. We observed the solar-mass Class protostar IRAS 16293-2422 with the heterodyne instrument HIFI in Bands la (480 - 560 GHz), 3b 

, ^, (858 - 961 GHz), and 4a (949 - 1061 GHz) as part of the Herschel key programme CHESS (Chemical HErschel Surveys of Star forming regions). 

Results. The deuterated form of the imidogen radical ND was detected and securely identified with 2 hyperfine component groups of its fundamen- 

I tal transition (N=0 -1) at 522.1 and 546.2 GHz, in absorption against the continuum background emitted from the nascent protostar. The 3 groups 

►^ I of hyperfine components of its hydrogenated counterpart NH were also detected in absorption. The absorption arises from the cold envelope, where 

I ' many deuterated species have been shown to be abundant. The estimated column densities are ~ 2 x 10'"* cm"- for NH and ~ 1.3 x lO''' cm"^ for 

^-f,. ND. We derive a very high deuterium fractionation with an [ND]/[NH] ratio of between 30 and 100%. 

V(»«, Conclusions. The deuterium fractionation of imidogen is of the same order of magnitude as that in other molecules, which suggests that an effi- 
_^ cient deuterium fractionation mechanism is at play. We discuss two possible formation pathways for ND, by means of either the reaction of N^ 
• with HD, or deuteron/proton exchange with NH. 

("^ Key words. ISM: abundances — ISM: molecules — stars: formation — stars: individual: IRAS 16293-2422 

o: 

^ [ 1. Introduction increase in the [H2D^]/[H3] ratio (and similar ratios involving 

• ^H multiply deuterated counterparts) enhances both the molecular 

^ The envelopes of low-mass Class protostars have been known deuterium fractionation by means of gas-phase reactions and the 

5^ to be characterised by high levels of molecular d euteration, a s atomic D/H ratio, therefore also increasing the de uteratio n of 

. .q , shown e.g., by Loinard et al. (12002]l, Roberts et al. (12007J|2002]l, species forming on grain surfaces (e.g. Roberts et al. l2U03T l. 
and the detection of multiply deuterated species (Ceccarelli et 

al. TOWfor D2CO, Parise et al. |2002| for CH D2OH, Parise et al. Mogj deuterated molecules detected to date are either satu- 
[2004 for CD3OH, van der Tak et al.|2002|for ND3). The high rated molecules or ions. The deuteration of radicals, which are 
abundance of deuterated molecules is believed to originate in important intermediate pieces leading to the formation of sat- 
the pre-stellar phase (Bacmann et al. ,2003): at the low temper- uj-ated molecules, has so far not been thoroughly investigated, 
atures (< 20 K) and high densities (n > 10^ cm^) prevailing in xheir observation could however provide an important clue to 
pre-stellar cores (and similarly in protostellar envelopes), heavy jhe formation process of the major deuterated species. For ex- 
molecules such as CO are depleted on the dust grains, so that ample, it is not yet completely understood if saturated deuter- 
species such as H2D+, D2H+, or DJ reach high abundances. The ated molecules are formed direcdy (from deuterated ion inter- 
mediates, whose dissociative recombination may lead to deuter- 



* Herschel is an ESA space observatory with science instruments ated radicals) like their hydrogenated counterparts, or if they are 
provided by European-led principal Investigator consortia and with im- formed by proton-deuteron exchange starting from the saturated 
portant participation from NASA hydrogenated molecule. The observation of the D/H ratio in the 
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radical may help us to disentangle these two processes. With this 
in mind, we study the deuteration ratio in the NH radical, which 
is formed from the dissociative recombination of NHj and NHJ, 
two precursors of ammonia in the gas-phase. Roueff et al. j2005l l 
argued that highly-deuterated ammonia can form from pure gas 
phase reactions. These models however require highly deuter- 
ated radicals, and we propose to observationally test this model. 

Light hydrides such as NH have fundamental transitions in 
the submillimetre domain: NH possesses 3 groups of hyperfine 
transitions around 950GHz-lTHz and its deuterated counter- 
part ND has its 3 groups of fundamental hyperfine structure tran- 
sitions around 490-550 GHz. At all of these frequencies apart 
from 490 GHz, the atmosphere is opaque and precludes ob- 
servations from the ground. However, these transitions can all 
be observed with the HIFI instrument on board the Herschel 
Space Observatory. Since hydrides generally have high Einstein 
Aui coefficients, only the ground state is expected to be signifi- 
cantly populated. As a consequence, it is possible to detect these 
species in absorption against the continuum background origi- 
nating from the heating of the nascent protostar. Moreover, as 
already mentioned above, deuterated molecules are very abun- 
dant in cold, dense medium, so that the envelopes of young pro- 
tostars represent ideal targets to study the deuteration of the NH 
radical. 

In the course of the CHESS key programme (Ceccarelli et al. 
1201 Oi l, we carried out a spectral survey of the frequencies cov- 
ered by HIFI in the young low-mass protostar IRAS 16293-2422 
(hereafter IRAS 16293). The frequency bands contain the funda- 
mental hyperfine transitions of ND and NH. In this Letter, we 
report the first detection by HIFI of ND, as well as the detection 
of the hyperfine transitions of NH. Section 2 describes the obser- 
vations and the determination of the radicals' column densities, 
and Sect. 3 discusses the implications of the derived deuterium 
fractionation on our understanding of deuteration processes in 
N-bearing species. 

2. Observations and results 

The solar-mass protostar IRAS 16293 was observed with the 
HIFI instrument (de Grauuw et al. 120101 Roelfsema et al . 1201 Oi l 
on board the Herschel Space Observatory (Pilbratt et al.. 120101) . 
as part of the HIFI guaranteed time key programme CHESS 
(Ceccarelli et al. l2010l l. A full spectral coverage of bands la (480 
- 560GHz), 3b (858 - 961 GHz), and 4a (949 - 1061 GHz) was 
performed on 2010 March 1, 19, and 3, respectively, using the 
HIFI spectral scan double beam switch (DBS) mode with opti- 
misation of the continuum. In this mode, the HIFI acousto-optic 
wide band spectrometer (WBS) was used, providing a spectral 
resolution of 1.1 MHz (-0.6 km s ' at 520 GHz and 0.3 km s ' 
at 1 THz) over an instantaneous bandwidth of 4x1 GHz. The tar- 
geted coordinates were a200() = 16'' 32'" 22?J5, (52ooo - - 24° 28' 
34.2". The DBS reference positions were situated approximately 
3' east and west of the source. The beam sizes at the frequencies 
of the ND and NH transitions are about 41" and 22", respec- 
tively, and the theoretical main beam and forward efficiencies 
are about 0.72 and 0.96, respectively. Calibration uncertainties 
are < 16 % for band 1 and < 32 % for bands 3 and 4 (June 2010 
estimates). 

The data were processed using the standard HIFI pipeline up 
to frequency and amplitude calibrations (level 2) with the ESA- 
supported package HIPE 2.8 (Ott et al. 12010). A single local 
oscillator tuning spectrum consists, for each polarisation, of 4 
sub-bands of ~ 1 GHz for the SIS bands (1 to 5). The 1 GHz 
scans are then exported as FITS files into CLASS/GILDAS for- 



may prior to data reduction and analysis using generic spectral 
survey tools developed in CLASS by our group. Spurious fea- 
tures in the spectra were first removed in each 1 GHz scan, and a 
low order (typically 3) polynomial baseline was then fitted over 
line-free regions to correct for residual bandpass effects. These 
polynomials were subtracted and used to determine an accurate 
continuum level by calculating their medians. Higher order os- 
cillations in the baseline are possible, but their amplitudes re- 
main low in the spectra considered here. The continuum values 
obtained in this way are well fitted by a degree 1 polynomial 
over the frequency range ~500-1200GHz. Sideband deconvolu- 
tion is computed with the minimisation algorithm of Comito & 
Schilke (2002) implemented into CLASS90 using the baseline- 
subtracted spectra. The single side-band continuum derived from 
the polynomial fit at the considered frequency (Table[T]) was 
eventually added to the spectra. The temperatures were con- 
verted to the Tnib scale, using the theoretical values of the main 
beam and forward efficiencies given above. 

We detect for the first time the 3 groups of hyperfine tran- 
sitions of ND at 491.9, 522.1, and 546.2 GHz of the N= - 1 
transition as well as the 3 groups of hyperfine components of 
the NH fundamental transitions at 946.5, 974.5 GHz, and 1 THz 
(see Klaus et al. ll997l for an energy diagram of the N= 0-1 tran- 
sition of NH). The transitions are seen in absorption against the 
continuum from the protostar and the hyperfine structure is par- 
tially resolved. The transitions of ND at 522.1 and at 546.2 GHz 
are presented in Fig.[T] From Fig.[Tl the lines can be unambigu- 
ously attributed to ND (but see Olofsson et al. 120071 who tenta- 
tively assign an unidentified line to ND). The transition of ND 
at 491.9 GHz is not shown here because some of its components 
are strongly blended with a bright H2CO line at 491.9683 GHz, 
with an SO2 line at 491.9347 GHz and with an HDCO line at 
491.9370 GHz. This transition was not used in our analysis and 
will not be discussed further. Figure|2] presents the 3 detected 
hyperfine groups of the fundamental transition for NH. 

The analysis was carried out by fitting the hyperfine struc- 
ture using the HFS method in the CLASS softwarqj. The line 
parameters (frequencies, Einstein Aui coefficients, and level de- 
generacies), originally from spectroscopic studies of Klaus et 
al. STWli and Saito & Goto (fT993b for NH and ND, respec- 
tively, were taken from the Cologne Database for Molecular 
Spectroscopy (Miiller et al. 120011 120051) . The routine assumes 
a zero-baseline spectrum (i.e., from which the continuum emis- 
sion has been subtracted) and therefore fits a line profile of the 
form (1 - e"'^)(7v(rex) - JviTbg) - T^), where r is the optical 
depth, Tc the intensity of the continuum, and Jv{Tex) and JviTbg) 
the radiation temperature for the excitation temperature Tex and 
the radiation temperature of the cosmological background, re- 
spectively. The parameters given by the HFS fitting procedure 
are reported in Table|2l Cols. (2) to (4). From the total optical 
depth T given by the fit, we can infer the opacities of single hy- 
perfine components by multiplying t by the relative intensities 
of these components. The total column density of the considered 
species is given by 



A^: 



StTV^ e(rex) 



gupAul 



e^ r 

hv I 

'tlcx — 1 J 



Tdv 



(1) 



where v is the frequency of the transition, c the velocity of light, 
^up the upper level degeneracy, Aui the spontaneous emission 



' http://www.iram.fr/IRAMFR/GILDAS 

- see the CLASS manual at |http://www.rram.es/IRAMES/otherDocuments/ma 
for more details 
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Table 1. Measured rms (at a frequency resolution of 1.1 MHz) 
and assumed continuum values for the different transitions (in 
Tmb scale). 
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Einstein coefficient, Tgx the excitation temperature, Q{T^y) the 
partition function at Tex, and E^^ the upper level energy. For the 
values of the excitation temperature and the frequencies consid- 
ered here, the factor with the exponentials is nearly equal to 1 . 
The velocity integrated optical depth of component / is given 

by TiAt) %L^ with t\ the optical depth of component / and Ay 
the linewidth determined by the HFS method (full width at half 
maximum). The obtained values of the column density are given 
for each group of hyperfine components in Col. (6) of Table|2] In 
our derivation of the excitation temperature, it was assumed that 
the zone where the absorption arises completely fills the emis- 
sion zone of the continuum. In the HFS fitting method, the de- 
termination of the optical depth comes only from the measured 
relative ratios of the hyperfine components and is independent 
of any assumption about the filling factor Overall, since the ex- 
ponential factor is close to 1 in EqlT] and the partition function 
depends only mildly on Tex (in the case of NH for example, it in- 
creases by 10% when Tex increases from 10 to 15 K), the inferred 
value of the column density depends very little on the assump- 
tions about the value of the continuum, beam efficiency values, 
or absolute calibration uncertainties. 

The column density found for NH is A?(NH) = (2.0 + 0.8) x 
lO'^^cm"^. For ND, the column density values found from the 2 
observed transitions differ slightly but are consistent within the 
error bars. A reasonable estimate for the ND column density is 
therefore: A^(ND) = (1.3 ± 0.8) X lO^^cmr^. 

With these values, we find a imidogen deuterium fraction of 
[ND]/[NH] ~ 30-100% in the envelope of IRAS 16293. We note 
however that because of the differences in beam sizes, the re- 
gion sampled by the ND observations is larger and on average 
characterised by lower densities and temperatures than the re- 
gion sampled by the NH observations for which the contribution 
of the warmer inner regions is larger. To take these effects into 
account, a more detailed radiative transfer modelling - taking 
into account the source structure and the exact coupling between 
the source and the instrument - would be necessary. The sim- 
ilarity of the determined linewidths (around 0.5 km s ' in both 
lines, when deconvolved from the resolution of the spectrome- 
ter) is consistent with the bulk of the signal from both molecules 
arising from the same (cold) region. 

3. Discussion and conclusions 

The very high deuterium fractionation observed in IRAS 16293 
is of the same order of magnitude as the molecular D/H ra- 
tio previously measured in other species in the envelope of 
IRAS16293, i.e., [HDCO/H2CO] = 15%, [NH2D]/[NH3] = 10% 
(van Dishoeck et al. 1995) [CH2DOH]/[CH3OH]=30% (Parise 
et al. I2004I I. The [ND]/[NH] ratio even appears to be the high- 
est measured deuterium fractionation in this source, though this 
result needs to be confirmed with a more detailed modelling. 
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Fig. 1 . ND (A^y = 1 - 1 1 ) transition at 546 . 1 5 GHz (upper panel) 
and ND {Nj = Oi - I2) transition at 522.1 GHz (lower panel). 
The filled histograms show the observed spectra and the solid 
line is the fit of the hyperfine structure performed using the 
CLASS software. At the bottom of each panel is a sketch of the 
positions and relative intensities of the hyperfine components. 
Components at the same frequency have been slightly put apart 
on the velocity axis for visualization reasons. The origin of the 
velocity axis was chosen so that the strongest component is at 
the LSR velocity of the source, v - 3.8kms '. The emission 
feature towards 522.12 GHz coincides with a methanol line. 



Our measured ratio is consistent with the gas-phase pre- 
dictions of Roueff et al. (I2005I I. although this model assumes 
that NH is a major product of the dissociative recombination of 
N2H +, a re sult that has since been refuted (Molek et al. 120071 
Talbi'20091). 

The main path for the gas-phase formation of NH in the 
dense interstellar medium is believed to be initiated by the en- 
dothermic (85 K) reaction N+ H- H2 — > NH+ + H (e.g., Galloway 
& Herbst [T989] l. which turns into NH+, NH+, and NH+ after 
successive additions of H2. We then expect NH to be mostly a 
product of the dissociative recombination of electrons with NHJ 
and NH3 (Thomas et al. 12005b . while NH^ leads to ammonia. 
One possible way of forming ND is by means of the reaction N+ 
+ HD — > ND^ + H. This reaction is less endothermic than the 
corresponding reaction with H2 ( 16 K- Marquette etal. 1988) so 
is more likely to occur at low temperatures. Furthermore, exper- 
iments indicate that the dissociative recombination of deuterated 
ions tends to eject H atoms preferentially (see Roueff et al. 120051 
and references therein). The dissociative recombination of e.g. 
NHD^ or NH2D^ should therefore favour the production of ND 
with respect to NH. On the other hand, the low abundance of HD 
(~ 3.2 X 10"^ [H2], Linskv l2007l l might make it diflicult to attain 
[ND]/[NH] ratios around our measured value. 

An alternative means of forming ND is by proton-deuteron 
reactions with NH, i.e., NH + H2D+ -^ NHD+ + H2, followed 
by a dissociative recombination of NHD^, which leads predom- 
inantly to the deuterated molecule. This scheme is similar to that 
proposed by Rodgers & Charnley (2001) to account for the high 
abundances of doubly deuterated ammonia, and relies on there 
being high abundances of species such as H2D^ . In cold dense 
regions depleted of heavy elements such as CO, the ratio of 
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Table 2. Parameters given by the HFS fit of each ND and NH transition, t is the total optical depth, Tex the excitation temperature, 
JviTex) the radiation temperature at the frequency of the transition, JviThg) the radiation temperature of the cosmological background 
at Tbg - 2.7 K, Tc the brightness temperature of the continuum, Au the linewidth (full width at half maximum), and A^ the total 
measured column density in the given species. Tex was deduced from Col. (2). Because of uncertainties in some stages of the data 
reduction process, we estimate that the linewidths are not known to better than about 0.25 MHz. It is this uncertainty that we quote 
here in Col. (3). 
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Fig. 2. NH (Nj = Oi - lo) transition at 946.5 GHz (top panel), 
NH (Ay = Oi - I2) transition at 974.5 GHz (middle panel) and 
NH (Ay = Oi - li) transition at lOOO.OGHz (bottom panel). The 
solid line fitting the spectra and the sketch at the bottom of each 
panel have the same signification as in Fig.[T] Note the HCN line 
towards 974.48 GHz. 



[H2D+]/[H+ ] is ind eed highly enhanced (e.g. Caselli et al. l2008l 
Caselli et al.'2003' Roberts et al.'2003^. The ion H2D+ has indeed 
been detected in the cold envelope of IRAS 16293 by Stark et al. 
(I2004I I. who estimated its abundance to be 2 x 10"'. Although 
both routes of ND formation considered here seem plausible, it 
is not possible for us to conclude anything about their relevance 
without a more complete chemical model. 

We have reported the first detection of the radical ND by 
the high resolution heterodyne instrument HIFI on board the 
Herschel Space observatory, towards the young solar-mass pro- 



tostar IRAS 16293. The hyperfine structure is seen in absorption 
against the continuum background from the protostar, for both 
ND and NH, which implies that the absorption arises from the 
cold envelope around the hot corino. The deuterium fractiona- 
tion measured is very high, with a ratio [ND]/[NH] between 30 
and 100%. Detailed radiative transfer and chemical modelling 
are needed to constrain the deuteration and formation of NH and 
this will be the subject of a future study. 
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